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We prove, via an elementary variational method, Id and 2d localization within the band gaps of a periodic 
Schrodinger operator for any mostly negative or mostly positive defect potential, V, whose depth is not too great 
compared to the size of the gap. In a similar way, we also prove sufficient conditions for Id and 2d localization 
below the ground state of such an operator. Furthermore, we extend our results to Id and 2d localization in d 
dimensions; for example, by a linear or planar defect in a 3d crystal. For the case of D-fold degenerate band 
edges, we also give sufficient conditions for localization of up to D states. 



I. INTRODUCTION 

Localization by impurities in periodic potentials with spec- 
tral gaps (band gaps) is a central topic in solid-state physics 
and semiconductor devices, 1 and it has direct analogues for 
other propagating-wave systems, such as for photonic crys- 
tals in electromagnetism. 2 We prove that, in one and two di- 
mensions (Id and 2d), localized solutions must arise in the 
gaps of a periodic Schrodinger operator for any "mostly" neg- 
ative or mostly positive defect potential whose depth is not 
too great compared to the size of the gap. To our knowl- 
edge, this is the first rigorous sufficient condition of this sort 
in 2d, aside from informal arguments based on effective-mass 
free-electron models close to quadratic gap edges, 3 extending 
an earlier theorem for Id localization in gaps, 4 and is quite 
different from the many rigorous asymptotic gap-localization 
results in the limit of very strong defect potentials. In 
addition to localization in gaps, we also prove Id and 2d 
localization below the ground state of a periodic potential 
for any mostly negative defect potential, extending earlier 
known results for localization in vacuum for any mostly neg- 
ative potentiafSHIH localization in t he per iodic case but for a 
strictly non-positive defect potentialji 4 ^ and Id localization 
for mostly negative defects. 4 Furthermore, we extend our re- 
sults to Id and 2d localization in d dimensions; for example, 
establishing localization for a linear or planar defect in a 3d 
crystal. For the case of £>-fold degenerate band edges, we 
show localization of D bound states for definite-sign defect 
potentials, and more generally relate the number of localized 
modes to the signs of the eigenvalues of a D x D matrix. Our 
proofs rely only on elementary variational eigenvalue bounds, 
generalizing an approach developed in Ref.[T3l 

One-dimensional localization in vacuum by an arbitrary at- 
tractive potential is easy to prove by a variational method, at 



the level of an undergraduate homework problem, 16 while a 
minimum non-zero depth for a potential well is required for 
localization in vacuum in three dimensions!^! The 2d case is 
more challenging to analyze. The first proof of 2d localiza- 
tion in vacuum by an arbitrary attractive potential was pre- 
sented in 1976 by Simon, 11 using more sophisticated tech- 
niques that also lead to asymptotic properties of the bound 
states. An el emen tary variational proof in vacuum was pro- 
posed by Pic q 1 18 ! 19 ! (and was adapted to Maxwell's equations 
in optical fibers by Bamberger and BonnePS), while a differ- 
ent variational proof was independently developed by Yang 
and de LlanoP^ An informal asymptotic argument utilizing 
properties of the vacuum Green's function was presented by 
EconomouPS On the other hand, in the case of periodic po- 
tentials, rigorous results for the existence of bound states from 
weak defect potentials V are more limited. Frank et al.^ an- 
alyzed the case of potentials V < localizing at energies 
below the ground state of a 2d periodic Schrodinger opera- 
tor; they not only proved that a localized state exists, but also 
bounded its energy in terms of a related vacuum localization 
problem. Here, we use a different technique, based on the 
variational method of Ref. 13 to prove existence of localized 
modes in Id and 2d for any indefinite-sign but "mostly" nega- 
tive defect potential V. This is closely related to our general- 
ization of Ref.[T3lto index guiding in the periodic Maxwell's 
equations! 2 ^ (Indefinite-sign localization was also considered 
for discrete Schrodinger operators. 21 ) As for localization in 
arbitrary gaps, however, we are not aware of published rigor- 
ous results in 2d that are valid for weak V. Prodan 4 proved 
that arbitrarily weak defects localize states in Id gaps, us- 
ing an asymptotic Birman-Schwinger technique similar to 
Ref.QTl imposing a mostly negative/positive condition on the 
defect potential identical to the one we use below; Prodan's 
result also applies to localization below the ground state of 
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a periodic Id potential. Various authors have shown that for 
strong defect potentials, those of the form XV where A^co, 
there exists a bound on the number of eigenvalues crossing 
the gap. 5 9 Localization has also been proved in the limit of 
high-order gaps in Id. 22 Another common, albeit somewhat 
informal, approach to gap localization is to consider local- 
ization for energies close to a non-degenerate quadratic band 
edge, making an effective-mass approximation and then quot- 
ing the results for vacuum. 3 Our proof of localization in gaps 
is non-asymptotic, does not assume a particular form of the 
band edge, and is an extension of the elementary variational 
technique of Ref. [13] for localization in vacuum. The trick 
behind the variational proof is to take the mid-gap energy, 
Eg, of the perturbed Schrodinger operator, H, and transform 
the question of energies near E g into an extremal eigenvalue 
problem. There are two typical ways to do this. One is to con- 
sider (H — Eg)^ 1 , which seems closely related to the Green's 
function method of Ref. [12] and to the Birmann-Schwinger 
condition, 23 but such an operator is hard to evaluate explic- 
itly and mainly lends itself to asymptotic analyses. Another 
method is to consider (H — Eg) 2 , recently used for another 
variational-type localization proof by Ref. [24] and it is this 
method that we adopt here. The same techniques are used in 
numerical methods to iteratively compute eigenvalues in the 
interior of the spectrum of a large matrix, where (H — Eg) -1 
corresponds to well-known shift-and-invert methods,^ and 
where (H — Eg) 2 has also been used (but is computation- 
ally suboptimal because it squar es the condition number) pSI2ZI 
Other possible technique s 9 ! 22 ! 28 ! have been suggested to us^ 
for proving such a theorem, but we are not aware of any 
published results for this problem other than the Id result 
of Ref. |4] Localization by weak defects is also related to 
self-focusing of solitons by nonlinearities, which was recently 
considered for spectral gaps in periodic potentials.^ 



II. BOUND STATES BELOW THE GROUND STATE OF A 
PERIODIC POTENTIAL 

Notation: In this section, unless otherwise stated, the sym- 
bol J f will stand for the integral over all of M 2 of the real- 
valued function /. 

The proof of the following will be a little simpler than the 
one in Sec.|TVl However, the ideas used here are almost the 



same as those of the proof in Sec. IV so it is hoped that af- 



ter going through this proof, the reader will easily follow the 
latter. Note that the theorem in this section has also been 
proved by other meth ods, b ut only for defect potentials that 
are strictly nonpositiveP^^This theorem will be slightly gen- 
eralized to allow for a defect potential that is localized in two 
dimensions, but has an arbitrary periodicity in all other dimen- 
sions, in Sec.|TTIl 



A. Problem statement 



Suppose we start with an unperturbed Hamiltonian, 



H Q = -V 2 + Vo, 



(1) 



where Vo is a periodic potential (possible generalization to 
non-periodic potentials is discussed in Sec. |V A) , which has 
a minimum energy eigenvalue Eq with at least one (degen- 



eracy will be explored in Sec. VB i corresponding "general- 
ized" eigenfunction ipo of the Bloch form: 1 a periodic function 
with the same periodicity as Vo multiplied by e lk x , giving a 
bounded \ipo\- We introduce a localized, indefinite (varying 
sign) defect, V, giving a new Hamiltonian H = Ho + V, sat- 
isfying the two conditions 



V\ip \ 2 < and 



(2) 



The contents of the individual sections are summarized as 
follows. Sec.|ll]provides a simple variational proof of the fact 
that any arbitrary mostly negative and sufficiently localized 
defect induces at least one bound state below the infimum of 
the spectrum of a 2d periodic Schrodinger operator. The Id 



lim V{r, 0) = 0. 



(3) 



case is proved in a similar fashion. Sec. Ill gives a generaliza- 
tion of that result by allowing the unperturbed Hamiltonian to 
be periodic in d dimensions and our defect potential to local- 
ize in two (or one) dimensions, but be periodic (with the same 
periodicity as the unperturbed potential) in the other d — 2 
(or d — 1) dimensions. Sec. |IV| presents our main results: it 
gives a sufficient condition for the existence of a bound state 
in a spectral gap of a 2d periodic Schrodinger operator. The 



results from Sec. HI can be applied to this case as well to al- 
low periodicity in any arbitrary extra number of dimensions 
(although the solution is only localized in two of these dimen- 
sions). The case of bound states in a gap confined along one 
dimension is proved in a similar way. Sec. [V] describes the 
case of degenerate band edges as well as some other possible 



generalizations. Sec. VI closes with some concluding remarks 
and possible future directions. 



Such a Vo and V are shown schematically in Fig.[T]along with 
a typical spectrum showing Eq (and also possibly a gap, which 
is considered in Sec. |IV) . We prove that conditions ([2]) and 
Q are sufficient to guarantee the existence of a bound state 
with energy lower than Eq. [Strictly speaking, in order for 
an energy < Eq to be a "bound" state, the potential V must 
be sufficiently localized so as to not alter the essential spec- 
trum (energies of non-localized modes) of Ho; ^ does not 
appear to be strong enough. Weyl's theorem provides a suffi- 
cient condition of relatively compact V, such as the Kato class 
of mostly square-integrable potentials, which seems to include 
most realistic localized defects. 3 ' 132 ] The proof is essentially a 
generalization of that in Ref. [131 



B. Proof of 2d localization 

Recall that the variational principle (or min-max princi- 
ple) states that the Rayleigh quotient (if;, Hip) / (ip, ip), where 
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FIG. 1: (Color online) Example of a periodic potential Vo and defect 
potential V with the spectrum a(Ho) located at the right. The spec- 
trum of the unperturbed Schrodinger operator has an infimum Eq, 
and may also have a gap in the set (Ei , E%), 



(-0, (j>) = J ip*(f> denotes the inner product between ijj and <f> 
(and ||t/>|| 2 denotes (ip,ip)), for any trial function ip (in the 
appropriate Sobolev space) is an upper bound for the ground- 
state eigenvalue of a Hermitian operator H ! 33 " 35 l Therefore, to 
prove the existence of an eigenvalue less than Eq (and thus a 
bound state), it suffices to find a ip such that 



Em 



(<P, (H - goW 



< 0. 



(4) 



The key point is to find a trial function that will work even for 
an arbitrarily small defect V. Motivated by Ref.[T3] we use the 
following trial function parametrized by a positive number a: 



i> — ^o7i where 7 



-(l+r)° 



(5) 



Once the appropriate trial function is selected, the remaining 
analysis is straightforward in principle — one simply plugs the 
trial function into E[ip] and finds an a > where it is nega- 
tive. The easiest way to do this is to take the a — > limit of 
the numerator: if this limit is negative, then there must also 
exist a small a > where it is negative. This process, which 
requires some care in taking the limits (limits and integration 
cannot be interchanged in general), is carried out as follows. 

Note that 7 is already in polar coordinates and that ip is 
normalizable for all such a since ipo is bounded. This trial 
function has the key physically motivated property that the 
limit of no localization, i.e. a — > 0, gives the unperturbed 
ground state V>o- We write down the first two derivatives of 7 
for future reference: 



7 



Or 



V7= -a(l + r) a ~ 1 r, 



(6) 



7" = |4 = a ( X + r )"~ 2 H 1 + r) a -a + l] 7. (7) 



When H — Eq acts only on ipQ m E[tp], the result is zero. 
The remaining terms in the Rayleigh-quotient numerator, de- 
noted by U[ip] , come from V and derivatives of 7. After some 
algebraic simplifications, U[ip] is given by (see the Appendix) 



v\4>\ 2 + I IV0I 2 



_V 2 ( 7 2 ) - 7 V 2 7 



Using equations |6]l and (FT), we obtain 



V 2 7 = a(l + r) c 



a(l + r) a - a- 



(8) 



(9) 



V 2 ( 7 2 ) = 2a(l + rf 



2a(l 



- r — a — 



7 



(10) 



Plugging these two formulas into U[ip] results in the concise 
form: 



U[ip] 



V\^+ I h/>o|^(l + r) 2Q - 2 7 2 



(11) 



Note that the denominator of E[i/j] (which is HV'II 2 ) is always 
positive and so does not affect the sign of U[ip]. We want 
to show that U[ijj], and thus E[tp], will be negative for some 
choice of a. This will be done by showing that the term on the 
right of equation (Hi tends to zero as a —> 0, while J V\ip\ 2 
will be negative inthis limit. Because | -0o 1 2 is bounded, we 
have 



|Vo|V(l + r-) 2Q - 2 7 2 



< 27rmax{|^ | } / a z {\ + r) Za - z Yr dr 



27rmax{|i/> | 2 }y aV Q " 2 7 2 (u - l)d% 



< 27rmax{|i/'o| 2 } 



7 du, (12) 



where we have made the substitution u = 1 + r and then 
bounded the integral again. Hence, it suffices to show that the 
latter integral tends to zero. We calculate this integral explic- 
itly via integration by parts: 



7 du 



-u e 



00 2 



,a-l„-2u c 



du 



3o 

4^' 



(13) 



Taking the limit as a —> yields zero as claimed above. 

The leftover J V\ip\ 2 term can be split into two parts. Let 
V = V + — V~ , where V + and V~ are the positive and neg- 
ative parts of V. Then we have 



(14) 



where each l^l 2 !/ 111 is a monotonically increasing function as 
a decreases. This allows us to use Lebesgue's monotone con- 
vergence theorem^ to interchange the limit with the integra- 
tion, arriving at 



lim / IVfV^ = / lim H 2 ^ = e~ 2 / iVol 2 ^ (15) 

a^aj J a^0 J 

for each part. Putting all the information together, we have 



Urn ( / V\^\ 2 + I \M'^{l + rY a - z r 



|Vo| 2 ^<0 (16) 



by our main assumption Hence, the variational principle 
says that there exists an eigenvalue below Eq for the system 
and so the theorem is proved. 
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C. Proof of Id localization 



The case for Id localization can be proved in an analo- 
gous way, with simpler calculations, by using the trial func- 
tion ipoe~ ax . 

A closely related result in Id was presented by Prodan: 4 
for any defect potential V satisfying and for any energy 
E < Eq, Prodan showed that there was some scaling XV 
with A > such that a bound state with energy E exists. 
Furthermore, the limit E — > E was shown to correspond to 
A — » 0, so that an arbitrarily weak potential satisfying Q must 
localize a state. 



to 



e "ctf 



• • • • %/e c ,9 



III. 2D LOCALIZATION IN d DIMENSIONS 

We would like to extend these results to an unperturbed 
Hamiltonian H Q = — V 2 + Vq that is periodic in d dimen- 
sions and where the defect potential V is localized along one 
or two dimensions and is periodic in the other dimensions. A 
physical example of 2d localization would be a linear defect 
or "quantum wire" in a 3d crystal, localizing a wavefunction 
to propagate along the line, whereas an example of Id local- 
ization in 3d would be a planar defect. 



A. Periodicity of Vb and V 

In this case, Vo is a periodic potential in d dimensions, 
while V is periodic only in d — 2 dimensions but localized 
in the other two. It is convenient to separate the periodic 
and nonperiodic coordinates of V by writing V as V(r, 9, z), 
where z £ ]R d ~ 2 , which is periodic in z with some lattice 
vectors. The first two coordinates are in polar form for con- 
venience in defining what we mean by "localization," which 
occurs in r only. Vo, on the other hand, is periodic in all di- 
mensions, but for convenience, we also write it as Vo(r, 9, z). 
A schematic 2d example of such a V and Vb, where V is local- 
ized in only one dimension and periodic in one other, is shown 
in Fig. [2] Note that it is irrelevant in our proof whether Vq is 
periodic in any (r, 9) plane, which only occurs if (r, 9) corre- 
sponds to a lattice plane 1 of Vb (such an orthogonal supercell 
exists only under certain conditions on the lattice vectors^. 

The corresponding Laplacian for these coordinates is: 

v 2 -v 2 +v 2 -^ + 1 ^ + l^+V^ 

V - V r , g + V z - I + + r2 I + 2^ g , • 

V 7 3=3 1 

(17) 

In this section and the next, the symbol J / will be inter- 
preted as 

f= rdr d9 d d - 2 zf(r,9,z), (18) 
Jo Jo Jn 

where il is the primitive cell in z. 



FIG. 2: (Color online) Linear defect for a periodic lattice (red cir- 
cles) in two dimensions is created by adding a defect potential (blue 
circles) that is periodic along the z direction, commensurate with the 
periodicity of the underlying lattice. This will localize a state with 
respect to the perpendicular directions, denoted r. 



B. Generalized localization below ground state 

Due to the periodicity of the potentials, we apply Bloch's 
theorem to reduce this problem to a unit cell in z, in which 
case a "localized" mode is a point eigenvalue as in the pre- 
vious section, and we can prove localization by an identical 
variational approach. 

In particular, all eigensolutions can be chosen in the Bloch 
form ipo = ipQ (r, 9, z)e tk ' z , where iJ)q is periodic in z with 
the same periodicity as V (and Vo) and k is some Bloch 
wavevector. 1 Substituting this into Hotp = Eijj, one obtains 
the "reduced" Hamiltonian 



k 2 + y , 



(19) 



whose eigenvalues E(k) are the k-dependent band structures. 
The domain is now only the primitive cell in z with periodic 
boundary conditions. For each k, there is a minimum en- 
ergy, Eq (k) , associated with each reduced Hamiltonian Hq , 
i.e. HqiPq = Eo(k)if;Q for some Bloch wavefunction ip*. 
(■0o elk z I s merely one of the Bloch wavefunctions for the un- 
derlying potential Vo such that the projection of the original 
c?-dimensional Bloch wavevector onto the z dimensions gives 
k. Therefore, i s bounded.) The claim we wish to prove 
is that the above conditions (just as in Sec. [II} guarantee the 
existence of a bound state with energy below E (k) for each 
k. Moreover, the previous theorem of Sec.|II]becomes merely 
the special case when k = and there are no extra dimensions 
of periodicity. 



C. Proof of 2d localization in d dimensions 



The proof follows the same outline as in Sec. [TT] (just re- 
place ip with ?/>q). There are only two small differences. Hq 
in ( 19 1 has the additional terms 2ik ■ V z and k 2 , but V z 7 = 
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so this additional term only acts on ip*, and is part of the 
(Ho — Eo)ipo = cancellation as before in <|8|> . The second 



difference is that these integrals are over I 



part over 



pd-2 



(where the 



is actually just over a 



over K. . This simply means that, in 



rimitive cell), instead of 
instead of factoring 



prirr 



out max{ | ipQ | 2 } in the bound, we factor out max{ | ipQ \ 2 } mul- 
tiplied by the volume of the primitive cell (since the remaining 
integrand is independent of z). 



D. Proof of Id localization 

The analysis is the same except we now have Vq = Vo(x, z) 
and ip = ip§e lk ' z e~ ax . 

IV. BOUND STATES WITHIN A BAND GAP 

Now, we suppose that there is a gap in the spectrum of Hq 
(at a particular k in d dimensions) from E\ up to E 2 (with 
corresponding band edge Bloch wavefunctions <ty\ and iptj) as 
shown in Fig. [T] and we prove sufficient conditions for V to 
localize a state along one or two dimensions with an energy in 
the gap (with periodicity as defined in the previous section). 
Intuitively, if V is mostly negative, then a state will be pulled 
down into the gap from the higher energy edge E 2 . If it is 
mostly positive instead, then a state will be pushed up into 
the gap from the lower energy edge E\. On the other hand, 
if the potential is too strong, then the band edge state will be 
pulled/pushed all the way across the gap and will cease to be 
localized, hinting that we may need some upper bound on the 
strength of the defect potential even if there is no lower bound. 
This will be made quantitative in Sec. |IV A] 

The idea behind the actual proof is to translate it to a sim- 
pler problem so that the methods of Sec.[ll]can be used. That 
is, we must somehow alter the operator H so that we are lo- 
calizing below the "ground state" of an altered H and can 
use the variational principle. This is achieved by considering 
(H - E g ) 2 , where E g = (E 2 - Ei)/2, instead of just H, 
which transforms the spectrum as shown in Fig. [3] This idea 
comes from the localization proofs of Re f. [24] as well as from 
numerical techniques .I2GI2ZI The trial function is motivated once 
again by Ref.fTSl 



A. Statement of the theorem 

For ease of notation, we will omit the k superscripts in this 
section. 

We will prove gap localization under the following condi- 
tions. The intuition that V should be mostly negative or posi- 
tive corresponds to our condition: 

J V\^pi\ 2 < and/or J V\tp 2 \ 2 > 0. (20) 

The intuition that V cannot be too big corresponds to: 



2\Ei — E„ 



> v 2 m\ 



(21) 



0(H) 



a(H-E g ) 



i 



E 2 Eg 



i 



(E -E g r 







E 1 Eg 



FIG. 3: The leftmost plot is a schematic picture of the spectrum 
o(H) (with some arbitrary origin 0). The first step is to shift the 
center of the gap to the zero position by subtracting E g . The mid- 
dle plot shows a(H — Eg). Then, the operator is squared and all 
energies become positive. The final plot shows a[(H — E g ) 2 ]: any 
eigenvalue introduced into the gap of the original operator will now 
be an extremal eigenvalue of the shifted and squared operator. 



where i depends on which case of (20i was satisfied; that is, 
from which band edge we are pulling a localized state. [Con- 
ditions pO])-([2T| ) can be merged into a single condition, be- 
low, that (28 i is negative.] We also require that V be suffi- 



ciently localized in r, corresponding to the condition: 



\V\\^ 



< oo 



and 



< oo. 



(22) 



These conditions are sufficient to guarantee a bound state in 
the gap for the perturbed Hamiltonian H = Hq + V, i.e. a 
localized state in the (r, 6) dimensions. 



B. Proof of 2d localization 

Considering one band edge at a time, we will prove local- 
ization of a bound state "pulled" from the band edge Ei for 
i = 1,2. The proof will be split up into several steps to make 
it easier to follow. The method employed is the same as that 
in SscM 

The variational principle in this problem can be used after 
shifting the center of the gap to the zero position and squar- 
ing the resulting operator thus making it sufficient to find a 
normalizable trial function ip such that 



&,(H-E g ) 2 TP) 



0M> 



<(E l ~E g ) 2 ^(AE) 2 



(23) 



or equivalently, 



\\(H-E g )ipf-(AEfM\\ 2 <0. (24) 

Consider the trial function ip = 7^, where 7 = exp{ — (1+ 
r) a } as before. Similar to Sec.|n] we will show that this ip will 
satisfy condition ( |24| for some small a providing conditions 
(|2"0"1)-(|2"2"1) are also satisfied. 

After some algebraic manipulation and using the fact that 
Hgipi = Eiipi, the left-hand side of ( 24 1 becomes (as shown 
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in more detail in the Appendix) 



\\(H - EM\- (^E) 2 H\\ 



(AE) J |^| 2 [V 2 (7 2 )- 2 7 V 2 7 ] 

+ 2(AE) J V\iP\ 2 + \\ViP\\ 2 

U[V|^| 2 -V( 7 2 )-2|^| 2 7 V 2 7 ] 
+ 4||V^-V 7 || 2 + ||^V 2 7 || 2 



2 J V 2 7 V|V 7 | 2 - V 7 . 



(25) 



(26) 



(27) 



It is a bit cumbersome, but we can separate the integrals into 
four main groups to avoid calculating each one individually. 



The first group consists of (25 i, but we have already shown 
that this integral tends to zero as a —> in Sec.[n] — see equa- 
tions <[12j» and ( fT3] >. 

The second group consists of all the integrals which involve 
V. These are all contained in (26 1. Taking the limit as a — > 



is allowed to pass through the integrals because everything is 
bounded by either a constant times V or a constant times V 2 , 
and our assumption ( |22| allows us to use Lebesgue's domi- 
nated convergence theorem 36 to commute the limits with inte- 
gration. The limit of these integrals as a — >• is 



V 2 \ip i \ 2 + 2{E i -E g )e 



(28) 



where the rightmost integrand of ( 26 1 vanishes for a — > be- 
cause differentiating 7 results in at least one a factor. Equa- 



tion ( |28] l is strictly negative under conditions (20i and (21 
Namely, the V 2 term is smaller than the V term by (21 
(20 1 implies that the (Ei — E g )V\ipi\ 2 integral is negative for 



and 



either i = 1 or i = 2. 

We now move on to the final terms: everything in ( |27"| i. We 
wish to show that they all tend to zero as a — > 0, and we 
can do this easily by concentrating on the term that decays the 
most slowly with respect to r, so that all the other terms which 
decay faster clearly go to zero as well (provided that the terms 
have the same or higher order for the a factor in front). The 
three integrals of (27 1, dropping bounded terms such as \4>i\ 2 , 
are [from left to right in ((27])] : 



(V 7 ) 2 = a 2 (l + r) 2 "- 2 7 2 
(V 2 7 ) 2 =a 2 (l + r) 2 "- 4 

V 2 7 V 7 = a 2 (l + r) 2 "- 3 



a(l + r) Q -a 
a(l + r) Q -a 



(29) 

2 

7 2 , (30) 
7 2 . (31) 



Each has at least an a 2 factor in front. Upon inspection, we 
find that the most slowly decaying term out of these three is 
the (V 7 ) 2 term, which goes as 1 jr 2 in the limit of a — > and 
its limit is 



lim 



as was already shown in ( fT~2| > and ( 13 i. Since (29 1 dominates 
p7j ), the other terms are all asymptotically bounded by some 
constant times this integral and hence must also have integral 
zero in the limit as a — > 0. 

What has been shown is that every term from the left-hand 



side of (24 1 vanishes except for the one term (28 i, which is 
negative. This establishes the existence of a bound state. 



C. Proof of Id localization 

The case for Id localization can be proved in an analo- 
gous way, with simpler calculations, by using the trial func- 

2 

tion TpQe~ ax just as before. 

A closely related result in Id was presented by ProdanPfor 
any defect potential V satisfying pDJ, and for any energy E in 
the gap, Prodan showed that there was some scaling XV with 
A > such that a bound state with energy E exists. Further- 
more, the limit as E approaches the band edge corresponding 

0, so that the limit 



a 2 (l + r) 2a - 2 -/ 2 r dr 



0, 



(32) 



to (20 1 was shown to correspond to A 
of an arbitrarily weak potential satisfying ( (20) ) must localize a 
state in the gap. 



V. SOME FURTHER GENERALIZATIONS 

A. Necessity of periodicity? 

In all our derivations, we did not actually explicitly use the 
fact that Vb was periodic in the dimensions where localization 
would take place. All we used were a few of the properties of 
periodic potentials. These are: 

• The energies are bounded from below. 

• There may be a finite gap inside the continuous spec- 
trum. 

• The generalized eigenfunctions corresponding to the 
infimum/gap-edge energies are bounded (and their 
derivatives are bounded). 

For gap localization, we also assumed that the squared 
operator was well-defined, making application of this the- 
orem to delta-function potentials (Kronig-Penney models) 
appear problematic, although it is possible that the diffi- 
culty is surmountable with a sufficiently careful limiting pro- 
cess. (In physical contexts, the difference between a theoret- 
ical infinite-depth Vq and a finite-depth approximation seems 
scarcely relevant.) Also, we assumed coinciding essential 
spectra for Ha and H in order to utilize the variational princi- 
ple. This means that there are some restrictions on how large 
of a perturbation V ca n be. However, Weyl's theorem states 
one sufficient condition) 3 ' l 32 l which appears to be satisfied for 
most physically interesting V's; see also the comments af- 
ter Q. 

The third assumption listed above may be more challeng- 
ing to prove for non-periodic potentials. We assumed that the 
energies at the infimum of the spectrum and/or the edges of 
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spectral gaps correspond to eigenvalues with bounded gener- 
alized eigenfunctions (^o or 4>i)- (This corresponds to the 
requirement of a regular ground state in Ref. [T5]) For peri- 
odic potentials, the existence of a band-edge solution for Vb 
follows from the well-known continuity of the band diagram 
as a function of the Bloch wavevector k. For non-periodic 
potentials Vq, however, this is not necessarily true. For ex- 
ample, for the Id half-well potential Vq(x) — oo for x < 
and = for x > 0, the eigenfunctions s'm(Kx) do not have 
a nonzero infimum-energy solution for n — > 0, and corre- 
spondingly it is well known that any perturbation V must ex- 
ceed some threshold depth before a bound state appears in that 
case — this is mathematically equivalent to the appearance of 
an odd bound-state solution ip{—x) — —ip(x) for Vq = and 
an even perturbation V(— x ) — V(x), which requires a thresh- 
old depth since the lowest-energy bound state in that case is 
even. It is not clear to us under what conditions the requisite 
infimum/gap-edge solutions exist for more general potentials, 
such as quasiperiodic potentials Vb, although some examples 
are given in Ref. [T5] 



B. Degeneracy at the band edges 

As mentioned earlier, it could happen that there are mul- 
tiple (degenerate) linearly independent ipi's corresponding to 
a given energy Ei at an edge of the gap and/or at the infi- 
mum of the spectrum. Our proof in the preceding sections is 
unaffected — there must still be at least one bound state local- 
ized by a suitable V, as long as the requisite conditions [(pi 



or (20 - 22 1] hold for at least one of the degenerate if>i wave- 
functions. Intuitively, however, one might expect to be able 
to prove a stronger theorem in this case — if Ei is D-fold de- 
generate, can one show that D localized states are formed by 
a suitable VI 

To prove the existence of more than one localized state, we 
can employ a generalization of the min-max theorem. For 
a single localized state, our proof involved the fact that the 
ground-state eigenvalue of a Hermitian operator O is bounded 
above by the Rayleigh quotient Q[ip] = (ip,(hp)/{ip,ip) for 
any ip. The generalization of this fact is that the n-th eigen- 
value A„ is bounded above byS^l 



A, 



< sup Q[ip], 
V>eS n 



(33) 



where S n is any n-dimensional subspace of the Sobolev space 
for O. We then wish to show that A„ < b for some bound b: 
O = H — Eq and 6 = for localization below the infimum 
of the spectrum, or O = (H - E g ) 2 and b = (E t - E g ) 2 
for localization in the gap from edge i. This is equivalent to 
proving that the Hermitian form B[ip, <p] = (ip, (O — b)<f>) 
is negative-definite for some n -dimensional subspace S n (i.e., 
B[V>,V] < for all -ip e S n ). 

If Ei is £>-fold degenerate, with degenerate generalized 
eigenfunctions ipj for £ = 1, . . . , D, then the analogue of our 
previous approach is to form the trial functions %jr = jtpf 
(whose span is a subspace Sd), compute the D x D matrix 
Bu> = B[ip e ,ip e ], and check whether it is negative definite 



as a — > 0. We wish to find the largest subspace S n of Sn for 
which B is negative-definite, which corresponds to the num- 
ber n of negative eigenvalues of B: this will be the number n 
of localized states that are guaranteed by the theorem. 

For localization below the infimum of the spectrum by a V 
satisfying ([3]), following exactly the same steps as in Sec. Ill 
proving that B[ip,<p] is negative-definite in this subspace re- 
duces to a generalization of condition Q. Specifically, show- 
ing B[ip,tp] < in the subspace for a — > reduces, via 
(16 1, to showing that J V\ip\ 2 < for every ip in the sub- 
space. In other words, the Hermitian form A[i/j, <j>] — (ip, V<j>) 
must be negative-definite in S„ . In the a — >• limit, this cor- 
responds to checking the eigenvalues of the D x D matrix 
Aw = A[ipg, ipQ ]: the number of negative eigenvalues of A 
is precisely the dimension of the largest negative-definite sub- 
space S n , and hence is the number of bound states that are 
guaranteed to be localized below the ground state of Vq. If we 
happen to have a strictly non-positive V < 0, then the Her- 
mitian form A is automatically negative-definite and we are 
guaranteed D localized modes. 

For localization in a gap by a V satisfying (|22]), pulling 
states from band-edge i, following exactly the same steps as 
in Sec. IV one finds that B[ip,4>] being negative-definite re- 
duces to a generalization of condition ( |28| : the Hermitian 
form G[ip, <j>] = {i>, [V 2 + 2{E l - E g )V}(f>} must be negative- 
definite in S n . As above, this simplifies for a — > to count- 
ing the number of negative eigenvalues of the D x D matrix 
Qui = G[ipf, ipf ]■ The number of negative eigenvalues is the 
number of solutions that are guaranteed to be localized from 
band-edge i. If V has sign everywhere opposite to Ei — E g 
and is sufficiently small (to overwhelm the V 2 term), then D 
eigenstates will be localized from this band edge. 

This analysis appears to be closely related to the asymptotic 
technique of Ref. |28l which also relates a number of bound 
modes to the number of eigenvalues of a given sign of a small 
matrix, via the Birman-Schwinger principle in the limit of 
weak perturbations, but that work only explicitly considered 
localization below the ground state of translation-invariant el- 
liptic and Schrtidinger-type unperturbed operators. 



VI. CONCLUDING REMARKS 

Although the existence of localized solutions from defects 
in periodic potentials and the effective-mass analogy with the 
vacuum case are well known as a practical computational and 
experimental matter, it is gratifying to have a general, explicit 
proof that localization in one and two dimensions occurs in 
a similar manner to localization in vacuum. A number of 
directions suggest themselves for future research. Although 
the simplicity of an elementary proof based on the min- 
max/variational theorem has its own appeal, the application 
of more sophisticated methods such as those of Ref. [TTI may 
reveal additional information about the nature of the localized 
state (such as its asymptotic localization length) that cannot 
be gleaned from a simple variational analysis. We would also 
like to transfer these results from the Schrodinger (scalar) pic- 
ture to the Maxwell (vector) one, in the context of localiza- 
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tion in band gaps of photonic crystals such as photonic-crystal 
fibers!^ A similar generalization to electromagnetism was al- 
ready obtained for localization below the infimum of the spec- 
trum (corresponding to total internal reflection in Maxwell's 
equations) for non-periodic 19 and periodic® media, and in 



photonic band gaps for sufficiently large defectsP^EH 

For localization in gaps, we should remark that the condi- 
tion ( f2T| i on the size of the perturbation V is somewhat un- 
satisfying. Intuitively, a "small" perturbation V could be one 
where either \V\ is small at every point or where \V\ is not 
small but the support of V is small. The latter case, however, 
of a large \V\ with a small integral, violates our smallness 



condition (21 1 for a sufficiently large \V\ no matter how small 
the support might be. This does not mean that there are no 
localized states in that limit — our proof only gives a sufficient 
condition for localization, not a necessary condition — but it 
suggests that some reformulation to handle this physically in- 
teresting possibility might be desirable. 
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Because U[tp] must always be real and 7 is assumed to be real, 
the imaginary part of tpgVtpo must integrate to zero. There- 
fore, we can replace 2ipQVip by tp^'VipQ + V'o^V'o m me 
integrand. We can also use the identity 7V7 = ^V(7 2 ) so 



that equation ( 35 1 becomes 



V\^\ 2 -\J V|Vo| 2 • V( 7 2 ) - J |V>o| 2 7V 2 7 . (36) 

We now rewrite middle term in another manner in order to 
eliminate the slowly decaying second derivative of 7, using 
integration by parts: 



V|^o| 2 -V( 7 2 )= / h/>olW)-y l^o| 2 V 2 ( 7 2 ) 

(37) 



l^o| 2 V 2 ( 7 2 ), 



(38) 



where " J g " stands for the boundary integral. The boundary 
integral is zero because the 7 term and it first derivative decay 
fast enough and tpo is bounded. Therefore, all we have to show 
in Sec.|n]is that 



VM 2 + ^ I |V>o| 2 V 2 ( 7 2 ) - / |^o| 2 7V 2 7 <0 (39) 



for some sufficiently small choice of the parameter a. 



B. Calculations for Sec.HVl 



Appendix 

Here we will provide more details for the calculations of 
Sees, [fl] and IV Let tp be an eigenstate multiplied by a real 
function 7 whose derivatives decay quickly enough, i.e. tp — 
?/>o7 where and 7 is real (we can think of the example used 
in this paper) and tpo is bounded and satisfies Hotpo = E^tpQ. 
We will explore some of the terms in both (tp, (H — Eo)tp) 
and ((H — E g )tp, (H — E g )tp). Our goals in the above sec- 
tions were to show that most terms tended to zero. In order for 
this to be apparent, we had to perform several algebraic ma- 
nipulations of the integrals that are given here more explicitly. 



A. Calculations for Sec.HTl 

We will derive the general formula used in Sec.|II]first. We 
have 

(H - E )tjj = (-V 2 Vo + V ijjah - E Q ^ ol 

+ Vtp- 2Vtp ■ V7 - VoV 2 7 
= Vip- 2V^o • V 7 - VoV 2 7 (34) 

since (-V 2 +y )V'o = E ^ , so that U[ip] = (tp, (H-E )ip) 



/ 7 VWo-V7- / \M"l^ Z l- (35) 



We will now derive the general formula used in Sec. IV 
The same assumptions hold as in the previous subsection, but 
with minor modifications from Sec.[lV] for example tp = tpcj 
where i = 1, 2 signifies the lower or upper edge of the band 
gap, respectively. First, we calculate (H — E g )tp. 

{H - E g )tp = (AE)tp + Vtp- 2Vtpi ■ V7 - ?/W 2 7, (40) 



where AE = E^ — E g as in Sec. IV Then the generalized 
equation for \\(H — E g )tp\\ 2 is given by (employing some of 
the trivial simplifications from the previous section) 



(A J B) 2 ||^|| 2 + 2(Aii;) / V\tp\ 2 + \\VtP\\ 



(41) 



fw\tpi\ 2 - V( 7 2 )-2 yV|^| 2 7V 2 7 (42) 
2(AE) J |^| 2 7 V 2 7 
||V^V 2 7 || 2 + 2 J V 2 7 V|^| 2 -V 7 
(AE) J V|^| 2 • V( 7 2 ) + 4||V^ • V 7 || 



(43) 
(44) 
(45) 



Notice that the first term in d45|l is similar to the middle term 



in ( 36 1 from the previous subsection, so the same analysis 
shows that (T45b becomes 



(46) 



+ (AE) / |^| 2 V 2 (7 2 ) + 4||V^-V 7 |^ 
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In ( |25)l o f Sec. |IV| this first term of (|46} is combined with 
term (|43|l. Several of the other terms are rearranged to make 



the presentation more concise. 
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